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. The lowest scalar and pseudoscalar glueball masses are evaluated by means of the time- 

dependent variational approach to the Yang-Mills gauge theory without fermions in the 
Ph ' Hamiltonian formalism within a Gaussian wavefunctional approximation. The glueball mass 

, is calculated as a pole of the propagator for a composite glueball field which consists of 

' ' ■ two massless gluons. The glueball propagator is here evaluated by using the linear response 

' theory for the composite external glueball field. As a result, a finite glueball mass is obtained 

I through the interaction between two massless gluons, in which the glueball mass depends on 

' ' , the QCD coupling constant g in the nonperturbative form. 

p., 

I , §1. Introduction 

CD ■ 

r"! ■ The hadronic and/or quark-gluonic world governed by mainly the strong in- 

teraction reveals very fruitful physics such as existence of various possible phases, 
' characteristic dynamical symmetries and various pattern of their symmetry breaking 

■ and so on in the hadronic and/or the quark-gluon matter.^) In the hadronic world 
which should be basically described by the quantum chromodynamics (QCD), the 

^Njj' ' color confinement is essential and the observed hadrons are mainly ordinary {3q)- 

■ baryons and (gg)-mesons. However, it is possible that the so-called color-singlet 
Tij- i exotic hadrons exist such as (3g(gg))-hadron like pentaquark hadron.^^ It is believed 

' that the color confinement occurs in QCD due to the non-Abelian nature^^ of the 

gauge interaction. This non-Abelian nature leads to the interaction between gauge 
fields themselves which represent gluons. Thus, it is interesting to consider a possible 
color-singlet state which only consists of gluons. This state is called the glueball. 
^ . The investigation about glueball has been performed widely in the theoretical 

side,^^ for example, by using the bag model, the flux tube model, QCD sum rule 
method, the lattice QCD simulation and so on. Especially, some lattice QCD simu- 
lations have given the glueball masses with certain spins and parties. On the other 
hand, in the experimental side,^^ some candidates of glueball states have been re- 
ported. However, the glueball states mix with the normal (gg)-meson states with 
the same spin and parity. Therefore, the definite information about glueballs such 
as masses and decay widths is not extracted experimentally until now. Thus, the in- 
formation about glueballs may be compared with the results obtained by the lattice 
QCD simulation. 

The many theoretical investigations about glueball have been carried out by the 
use of the effective model of QCD. In their several investigation, the finite gluon 
mass is assumed in a certain treatment and the background of the model is not 
so clear because the gluon does not appear explicitly in a certain model. Thus, it 
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may be necessary and interesting to deal with the gluebah starting from the QCD 
Lagrangian itself. 

In our previous papcr^^ which is rcfcrccd to as (I), the time-dependent variational 
method has been formulated for the Yang-Mills gauge theory without fermions in 
the Hamiltonian formalism, in which a Gaussian wavefunctional has been adopted 
as a possible trial state. This method presents an approximate treatment for the 
Yang-Mills gauge theory within the Gaussian wavefunctional approximation which 
corresponds to the Hartree approximation. Further, by the help of the linear re- 
sponse theory,'^) it may be possible that an approximation beyond the Hartree ap- 
proximation, which may correspond to the random phase approximation (RPA) in 
many-body physics, is obtained in a non-pertTirbative way. In Ref.8), it was shown 
that the Goldstonc theorem is satisfied in the time-dependent variational approach 
to the linear sigma model by the help of the linear response theory, while the Gold- 
stone theorem breaks down in the tree- level approximation. Further, in Ref.9), the 
pion and sigma meson masses have been calculated by using the linear response 
theory in the linear sigma model. The same approach is possible to calculate the 
glueball masses in the QCD without quarks. Since the time-dependent variational 
method may be suitable for the use of the linear response theory in the quantum field 
theory, the Hamiltonian formalism is applied. In this paper, starting from the QCD 
Lagrangian density without quarks, the lowest scalar (0+) and pseudoscalar (0~) 
glueball masses are investigated. Then, the glueball masses are obtained reasonably 
compared with the results of the lattice QCD simulation. The glueball in this pa- 
per consists of two massless gluons which interact by the self-interaction due to the 
characteristic feature of QCD. As a result, through the interaction, the glueball gets 
mass. 

This paper is organized as follows: In the next section, the time-dependent 
variational approach to the Yang- Mills gauge theory without fermions is summarized 
following (I). The formalism partially owes to Ref.lO). In §3, the method to calculate 
the glueball mass is explained. In §4, the lowest scalar and pseudoscalar glueball 
masses arc evaluated in the modified minimal subtraction scheme which is described 
in detail in Appendix A. Also, the dependence of the glueball masses on the QCD 
coupling constant is shown. The last section is devoted to a summary and concluding 
remarks. In Appendix B, it is shown that the decay width may appear if an imaginary 
part of a response function investigated in this paper is considered seriously. In 
Appendix C, it is verified that the gluon mass is zero under the approximation used 
here in this formalism. 

§2. Recapitulation of the time-dependent variational approach to the 
Yang-Mills gauge theory without fermions 

In this section, the time-dependent variational approach to the su{Nc) Yang- 
Mills gauge theory, which has been developed in our paper, is given and summarized 
following (I) to make this paper be self-contained. In order to formulate the time- 
dependent variational method, the Hamiltonian formalism is adopted,^^) in which a 
Gaussian wavefunctional is applied as one of possible trial states. This trial state 
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includes a mean field and quantum fluctuations around it as variational functions. 
The time-development of the mean field and quantum fiuctuations can be described 
under this Gaussian approximation. The Gaussian state used here corresponds to a 
squeezed state. 

2.1. Hamiltonian formalism of Yang-Mills gauge theory without fermions 

As is well known, the su{Nc) Yang-Mills gauge theory leads to a constrained 
system. Therefore, it is necessary to impose a constraint conditions. As was discussed 
in (I), as a result, the Hamiltonian density can be simply expressed^^^ so as to be in 
Eq.(L2-10). Namely 



{E''f + {B''f , (2-1) 



where roman letters a denotes color indices. Here, we define = {E^, Ey, E^) and 
so on and 

E"" = -i 



= V X A" - ^gf'^'A'' X A" , (2-2) 

where A"' is a gauge field and its conjugate field is identical with a color-electric field 
E"'. Also, the color-magnetic field is introduced and /"'"^ is a structure constant 
of su{Nc) Lie algebra. In the functional Schrodinger representation^^) for a field 

theory, the conjugate momentum field E"" in the gauge theory is represented as a 
functional derivative with respect to a gauge field A", in which they obey canonical 
commutation relations: [ A1(x), E^-{y) ] = i8ij5abSi^{x — y). 

We formulate the time-dependent variational method for the Yang-Mils gauge 
theory. Then, it is necessary to introduce a trial state \^) for variation. Here, in this 
paper, the trial wavefunctional is adopted as a Gaussian form as follows: 

#(A«) = (A-^l^) 

= Ar-^exp(z(:E|A- A))exp^-(A- A|^-zi:|A- A)^ , (2-3) 

where A/" is a normalization factor and we use abbreviated notations as 
(E\A) = j d^xE^ix^t) ■ A"(a;) , 

(^1^1^) = 11^'^ Anx)\G-'t^ix,y,t)A''^{y) . (2-4) 

Here, A'^{x,t), E^{x,t), Gfj{x,y,t) and Efj{x,y,t) are the variational functions 

in which A'^(x,t) and E'^{x,t) correspond to the expectation values of the field 
operators Af(x) and Ef{x), respectively, and two-point functions Gfj{x,y,t) and 
Efj{x,y,t) are related to the expectation values of the composite operators as fol- 
lows: 

mf{xm = A^{x,t), 
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{<P\Et{x)\<P) = Et{x,t) , 

= A':ix,t)A';iy,t) + Gt^ix,y,t) , 

{^\Ef{x)E^{ym = E'^{x,t)^{y,t) + ^G-'t^{x,y,t)+4{EGE)t^{x,y,t^ , 

(#|^f(a;)£;,^(y)|#) = A"(a;,t)E;(y,t) + 2(Gr)«f(a;,y,i) . (2-5) 

Thus, it is understood that ^(x, t) represents a mean field for the field A^{x) and a 
diagonal element of two-point function Gfj{x,y,t), namely Gf^{x,x,t), represents 
a quantum fluctuations around the mean field A'^{x,t). The state \$) in (2-3) is 
identical with the squeezed state. Thus, in this paper, the lowest excitation mode 
with a certain quantum number around vacuum is only treated.^) 

To determine the time-dependence of the state |<P) or a Gaussian wavefunctional 
it is necessary to determine the time-dependence of the variational functions 
A^{x,t), Ej^{x,t), Gfj{x,y,t) and Ufj{x,y,t). The time-development of the state 
under the Hamiltonian density H is governed by the time-dependent variational 
principle in general: 

6 J dt{^\i^ - J d^xn\<P) = . (2-6) 
Here, the Hamiltonian derived from the Hamiltonian density (2-1) can be expressed 



as 



{Hq) = j £x{$\no\^) 

= J d^x(^^Wix)-B\x) + ^E\x)-E\x) + ^Tt{x\G-^\x) 

+2Tr(a;|i;Gr|a;) + -Tr{x\KG\x) + ^ (TrlSiT" {x\G\x)]f 
2 8 

+^Tr [SiT''{x\G\x)SiT''{x\G\x)A , (2-7) 



where 



^ = ^ijkdjX - lar'^^jkA^Al , 
K = {-iS ■ Df - gS ■ B , 

D = V-igA, Ai=T,T\ Bi = B'^T''. (2-8) 



Here, eijk is a complete antisymmetric tensor and S implies a spin 1 matrix. The 
time-dependent variational principle (2-6) leads to the following equations of motion: 
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where the dot represents a time-derivative, that is, A {x, t) = dA {x,t)/dt and so 
on. These equations of motion are identical with the canonical equations of motion 
since the canonicity conditions, which are developed in the theory of collective 
motion in nuclei, for the variational functions are implicitly imposed. Under the 
Hamiltonian density Ho, a possible set of solutions with respect to A and E is 
given by 

A''{x,t) =E''{x,t) = . (2-10) 

Instead of the equations of motion for the two-point function Gfj{x,y,t) and 
Ufj{x,y,t) in Eq.(2-9b), we can reformulate the equations of motion by introducing 
the reduced density matrix. The reduced density matrix M. is defined as 



{Ef{x,t)E]{y,t)) i{Ef{xM]{y,t)) - \ ) 

-2iiGU)t^{x,y,t) G-^ix,y,t) 
I iG-^)t^ix, y, t) + 4{EGE)-^ix, y, t) 2i{SG)^^ix, y, t) 



,(2-11) 



where and Ef represent the quantum fluctuations around mean fields, which are 
defined by 

AUx,t) ^ A?{x) - (AUx)) = A?{x)-^{x,t) , 

Et{x,t) ^ Efix) - {Et{x)) = Etix) El{x,t) . (2-12) 

As for the reduced density matrix, it has been shown that the time-development of 
the reduced density matrix is governed by the following Liouville-von Neumann type 
equation of motion: 

iM'i^{x,y,t) = [no, M ]t^{x,y,t) , (2-13a) 
where the Hamiltonian matrix Ho corresponding to Ho is introduced as 

H^hx, y,t)=( ]s\x-y), (2.13b) 



rrnx,t) 



^2 



r^fix, t) = Kt^ + {SkT%x\G\x)SkT% + ^ {SkT% TV [S^T- {x\G\x)] . 

(2- 13c) 

The reduced density matrix (2-11) satisfies the following relation: 

M^=(^l i) ■ (2-14) 

Therefore, it is concluded that the eigenvalue of the reduced density matrix M. itself 
is -1-1/2 or —1/2. Let the eigenstate for M. be |(7l/2, n, a, i) where cr = it and 
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n represents a certain quantum number. Thus, the eigenvalue equation for M. is 
written as 

Mt^\al/2,n,b,j) = a^\al/2,n,a,i) . (2-15) 

Then, let us introduce the mode functions u and v in the coordinate representation 
by using the eigenstates for Ai: 

M+l/2,„,c.,i)=(<(jJj) . (2.16) 

Then, it is possible to express the reduced density matrix in terms of the above mode 
functions, namely, the spectral decomposition can be carried out: 



n(<T>0) 



+ (-W))<-""5«'""5« 



(2-17) 



where /„ = 1/2. Here, the mode functions satisfy the following orthonormalized 
conditions as 

^ J {v*^,f{x)Unt{x) + ul,f{x)Vnt{x)) = S^n' , 

a,i 

^ I Sx «,«(a;)<«(a=) - <,«(a.)<«(a;)) = , 

a,i 

j {Un't{x)Vnt{x) - Vn'tix)Unt{x)) = . (2-18) 

a,i 

For the equation of motion (2-13a), if the reduced density matrix Ai = Mq has 
no time-dependence, the equation of motion is reduced to [ "Hq ) -^0 ] = 0, namely, 
the reduced density matrix and the Hamiltonian matrix is commutable each other. 
In this situation, it is possible to diagonalize Mo and Ho simultaneously. Thus, the 
diagonal basis \aika) can be introduced, where the following relations are satisfied: 

Holaika) = E^\aika) , Mo\aika) = f^\aika) , (2-19) 

E^ = aEk, /k = ^ • ^ , ^ = ± • 

Here, the second equation in Eq.(2T9) is identical with Eq.(2-15), so the diagonal 
basis is nothing but |o"l/2, n = k,a,i): 

\aika) = \al/2,k,a,i) . (2-20) 

The Hamiltonian matrix and the reduced density matrix are diagonalized by using 
the mode functions in the momentum representation such as 



U-'noU=[ ^ _^^) , U-'MoU=\^f^ _i ) . (2-21) 
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Here, a unitary matrix U in the momentum space is obtained as 

_ ( -(^)?| C/-Vfc)«^-^*(^)f ^^%\ (2-22) 



where 



U{k)1^ = U^tSijSab , V{k)1^ = VkfSijSab , (2-23) 



(< ,<) = (fc|l/2,fe,a,0 
= |fc| . 



1 



(V2^ 




In the diagonal basis, the Hamiltonian matrix and the reduced density matrix are 
expressed in the forms of the diagonal matrices as (2-21). Originally, both arc ex- 
pressed in the form of Eqs.(2-13b) and (2-11), respectively. The relation between 
the original basis {fo;)}, which we implicitly used in Eqs.(2-13b) and (2-11), and the 
diagonal basis {\aika)} is given by the relation 

||a) = U-'^\aika) . (2-24) 

§3. Scalar and pseudoscalar glueball masses 

In this section, a scalar and pseudoscalar glueball masses are calculated in the 
framework of the time-dependent variational method within the Gaussian approx- 
imation developed in the previous section. In general, the propagator S can be 
obtained by using the generating function of the connected Green function, W[J], 
as 

where J represents a source current. Here, since the expectation value of the field 
operator (p is obtained as ip = dW[J]/SJ, the propagator can be expressed as 

Su^'^. (3.2, 

where the subscripts / and J imply certain indices. In our Hamiltonian formalism, 
we introduce the external field (pi{x) with a source current Ji{x), and the external 
Hamiltonian H^xt should be added to the Hamiltonian Hq: 

H^^t = J d^xJi{x, t)(pi{x) = ee-'^^ J d^xe^'^ ''Ji^i{x) , (3-3) 

where the source current Ji{x) is assumed as Ji{x) = ee~'*^J/. Thus, as for the 
expectation value of the field operator ipi which we write ((^7), the same time- and 
coordinate-dependence as that of the source current is expected: 

{ipi{x)) = Pie-"-*+'^-^ . (3-4) 
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Therefore, the propagator Sjj is evaluated by 

Sij = Urn 4- • (3-5) 

e^O ejj 

Namely, a mass of particle represented by the field operator ipi is calculated as a 
pole of the propagator, namely, the pole of limg^o 

Under the existence of the source current Jj, the reduced density matrix M. is 
shifted from Mq which means the reduced density matrix without the source term. 
Of course, the Hamiltonian matrix T-L is also shifted: 

n = no + sn , 

M=Mo + dM , (3-6) 

where the quantities with subscript represent those without the source term. Since 
the change of the reduced density matrix leads to the change of its (l,2)-component, 
G'^j in Eq.(2-ll), then the Hamiltonian matrix which contains G^j* is also changed 
through r^f in Eqs.(2-13b) and (2- 13c), namely F^f F^f + SF^f. Thus, up to 
the order of e, the Liouville von-Neumann type equation of motion for the reduced 
density matrix can be expressed as 

16^4 = 1^0 , SM] + [ (5Hext , A^o ] + [ SHind , Mo] , (3-7) 

ray - f ^ ^''i^"-^ 

Here, the shift of the Hamiltonian matrix, SH, is divided into two parts^ namely ST-L = 
Sy-ind + STiext, in which STiiad represents the shift due to SF and 5Tiext represents 
the shift occurring directly from the introduction of the external Hamiltonian -ffext- 

3.1. Scalar glueball mass: 0^ 

The scalar glueball field is constructed as the following composite operator: 

VI ^ \K^Fr = EtEt + BfBf . (3-8) 

As for the external field B^Bf, the response 5Ai has the same form as that for EfEf 
in our Gaussian approximation. Thus, we adopt the external Hamiltonian simply as 
follows: 

i^ext = Jd^^l d'yJt^{x,y,t)Et{x)E]{y) , (3-9) 

where the source current J^j{x, y,t) is proportional to SijSabS^{x — y) for 0"^ glueball 
state. We can derive the Liouville-von Neumann type equation of motion for the 
reduced density matrix as is similar to Eq.(2-13a). 

iMf^ix, y,t) = [Hj, M ^^{x, y, t) , (3-10) 
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where 

By using the equation [ Hq , A^o ] = and, first, omitting the induced term 
STiindj the foUowing equation of motion is obtained: 

iSM = ['Ho , SM] + [ SHe^t , Mo]. (3-12) 

We call SM derived by neglecting the induced term the bare response. 
Let us take the source current in the form 

J^j'ix, y, t) = j|^e-''^*+^'i-^53(^ _ _ ^3.;l3) 

Since STicxt is proportional to e~*'^*+*'?'^, then, the bare response 6Ai is also propor- 
tional to e~*'^*+*9'^, which leads to i6Ai = uj6Ai. Thus, by using the diagonal basis 
{\aika)} which obeys the eigenvalue equations in (2-19), Eq.(3-12) is recast into 

{aika\SM\bjk'a') = lf-^A^{aika\6ne^\bjk'a') . (3-14) 

Here, we can rewrite the matrix elements in terms of the original basis as 
{aika\6ne^t\bjk'a') = Y,{<^ika\a)) {{a\5'H^^\l5)) mhjk' a') 

'0 2 J,^ \ 



(A;|e-^«^|fe') 



mj 

(j3(fc'-fc + g)e-^'^*.(3-15) 



As is similar to the above transformation, inversely, the shift of the reduced density 
matrix SM. can be expressed in terms of the diagonal basis as 

( {elk+\5M\mfk'+) {elk+\6M\mfk' -)\ rj-i^,i-.fb 



^^'^^'^^y{elk-\6M\mfk'+) {elk-\5M\mfk' -) ) ^ )mj • (3-16) 



Thus, we obtain the shift of G^j and so on such as 



(5(i«4) = ui,f{aik+\SM\bjk'+)ul,^j + ul^{aik-\5M\bjk'+)ui,''j 
-Ukf{aik+\SM\bjk'-)u^>'' - ulf{aik-\6M\bjk'-)uk, 
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= 59+4^ {k, k - q)5\k' -k + q)e-'^* 
SiEfE^j) = Vkt{aik+\SM\bjk'+)vl,''j - v^'^{aik-\5M\bjk'+)vl,''j 
+Vk^{aik+\5M\bjk'-)vi,,^j - vl'^{aik-\6M\bjk'-)vk, 



= 6g-+^^{k, k - q)5\k' -k + q)e 



3nJ 



-iujt 



(3-17) 



with (3-14) and (3-15). 

Next, let us include the induced term SHind- First, the shift of Gfj{x,x) is 
calculated from Eq.(3-17) as 



SGf^{x,x) = {x\5{A^A'')\x)= I d'k I d''k'{x\k){k\S{AfA'')\k'){k'\x) 



\a Ab\ I r,/\ / I,/ 1 



= J (fk j (fk' 



„i(k-k')-x 

(27r)3 



{k\5{A<lA'm 



Sk 



/ 



Sg+.t^{k,k-q)e 



-iuit 



(3-18) 



Similarly, 

5Sf^{x,x) = {x\S{EfE^j)\x) = Jd^kJ d^k' {x\k){k\S{Ef&j)\k'){k'\x) 



J d^k J d^k' 



,i(k-k')-x 



(27r)2 



-{k\5{EtE\)\k') 



<fk 



j^6g.4^ik,k-q)e-^-' 



(3-19) 



Thus, the shift of F^J' in the Hamiltonian matrix, Sr^J" in Eq.(3-7), is given in the 



momentum representation as 

{k\5r^f\k') = 6r,f{oj, q)6^{k' -k + q)e-^^ 



d^k 



{SkT'-6g+4k,k-q)-SkT 



:c\ab 



+ ^-{SkT-)1^Tv [SkT'.6g+.{k,k-q)] 



. (3-20) 



Here, from Eq.(3-18), we obtain a^^{q) = J ^Sg+-^^{k,k- q). Thus, dE^fiq) 
can be expressed as follows: 



(3-21) 
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Further, the induced term of the Hamiltonian matrix is obtained in the diagonal 
basis as 



{aika\SHind\bjk'a') = U-\k)tt 



mj 



-u(k)-5r(g)f^n(k')/4 



3 "k'j 



5^{k' -k + q)e 



-iuit 



(3-22) 



Thus, from the equation of motion (3-7), we can obtain the response SM. in the 
diagonal basis in the same way that the bare response was derived in Eq.(3-16): 



{aika\SM\bjk'a' 

f<j' fa 



{aika\5Tiir,d\bjk'a') + {aika\5Tiext\bjk'a') 



k 

* a 



Vki 



* b 
"k' 



x5^{k' -k + q)e 



-iojt 



(3-23) 



Now, let us return to evaluating the scalar glueball mass. Since the field operator 
representing the scalar glueball is in Eq.(3-8), the source current (3-13) should be 
applied as follows: 



ab 



ij^ab ■ 



(3-24) 



In this case, the structure of the Lorentz and color indices of the response is the 
same one as J^j' such as afj oc SijSab- Namely, 



(3-25) 



Then, the shift SF^f' in Eq.(3-21) is simply written as 



Sr,f{u;,q) = 9\SkUTr'aiiSk)ij{n 



c\eb 



(3-26) 



Here, we used {Sk)ii = 0, {T'^Y^ = and J2^(t^T' {T'r = c^Sab with cg = 3 for 
s'u(3)-algebra. Further, it should be noticed that /+ = 1/2 and /_ = —1/2, which 
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leads to {aik+\6M\bjk'+) = {aik-\6M\bjk' -) = 0. Substituting (3-17) into (3-18) 
with {aika\5M.\bjk'a') obtained in Eq.(3-23), we can express in terms of SF^j' 

and J^j' or e. However, SF^J' is also expressed by aj in Eq.(3-26). As a result, we 
obtain 

miq) = 9^C3 ai{q) - ai{q)j flo{u, q) - eEoiuj, q) , (3-27) 



where we defined 



<fk 1 1 



ni^\u,q) 



d^fc 1 1 



Thus, we can solve Eq.(3-27) with respect to aj(g): 

= 4^"'^^ . (3.29) 
252c3iTo(a;,q)-l 

In the same way, we can obtain j3i{q) as 

^i^l) = 7 Wc^hio^^qf + ^o(a;,g)(l - 2^203770(0;, g))! ,(3-30) 

1 - 2g^C3no{uj, q) L J 

^0(0;, q) = S^"^ {uj, q) - E^'^^ (w, q) , 

d^k ^k^k-q 1 



St\u:,q)= J 



(27r)3 2 co±{Ek + Ey,_^) 



Thus, the scalar glueball propagator Sff'{co,q) is obtained in the general manner 
given in (3-5): 

5fi"(g) = lim^ . (3-31) 

From the relation (3-31) with (3-30), the pole of propagator Sf^{q) gives the scalar 
glueball mass Mq+: 

l-2g^C3no{Mo+,0) = . (3-32) 

3.2. Pseudoscalar glueball mass: 0~ 

The field operator of the pseudoscalar glueball is given by 

VI = \F;.Fr = \ {EtBt + BfEt) , (3-33) 
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where Fa'^ = (l/2)e^'^'^'^ Fp^. Thus, let us consider the external Hamiltonian with 
external current Jf-^ as 

iJext = d'yJ^J'{x,y,t) ■ \ {Et{x)B){y) + Bt{x)E]{y)) . (3-34) 

Here, the relation J^j{x,y,t) = Jj!^{y,x,t) is satisfied. As is similar to the case of 
the scalar glueball, the external Hamiltonian matrix in the equation of motion for 
the reduced density matrix is obtained like Eq.(3-ll): 



6n^^ttj{x,y,t) 





'jfjix^y^t) 



(3-35) 



where we define 



jf^{x,y,t)^{S)u'-Tjf^{x,y,t) , 



'iy\'"itfi"j- \-^)3'J ■ (3-36) 

Here, ^ (^^) means the derivative with respect to x (y) for the quantity on the 
right-hand (left-hand) side. As the same way deriving the external Hamiltonian 
matrix in the diagonal basis in Eq.(3-15), we obtain 



{aikalSHextlbjk'a') 



"ki 



x5\k' -k + q)e-'''^ , (3-37) 



where 



JtJ'{x,y,t) = J d\Jt^{q)e-^-'e'^-^5\x-y) . 



(3-38) 



Including the induced term dUiad-, the shift of the reduced density matrix is calcu- 
lated in the same manner deriving Eq.(3-22). As a result, we obtain 

{aika\5M\bjk'u') 



fa fa 



+ 



"ki 



x5%k' - k + q)e-''^^ . (3-39) 



As for the pseudoscalar glueball, let us take J^j' in the form 



4^ = eSi.S. 



(3-40) 
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Under this form, the source current appearing in Eq.(3-37) has the foUowing form: 

Jtfil) = ie{S)ij ■ qSab = -eekijQkSab , (3-41) 

The response a°^^{q)'m. Eq.(3-18) has the same structure with respect to the Lorentz 
and color indices as that of Jif{q), namely 

afjiq) = OL ■ {S)ijSab = akiq)eijkSab , (3-42) 

where ak{q) is introduced. Following the same procedure in which ai{q) in Eq.(3-27) 
was derived, we can obtain the following relation as 



c^^iq) = Sr,f{q)no{q) + (j^t^{q) - J,f{q)) To{q) , (3-43) 
where 



5r^Hq)=g^c^af'{q) (3-44) 



and we defined 



To{q)=T^-\q) + n-'\q), 
1 



Thus, ai{q) which was introduced in Eq.(3-42) is determined: 

2eqkfo{q) ar\ 
ah = . (3-46) 

l-g^CsHoiq) 

Instead of (3-19), we should investigate the following quantity: 

^<5g__t(fc,fc-q)e--* 
= ejimqie-"'''Pt{q) , (3-47) 

where 

P'l 



/d k r 
[vk^{aik+\6M\bmk - q+)u%_ 



b 

-qm 



-vli{aik-\SM\bmk - q+)ul_g^ - V]^f {aik+\SM\bmk - q-)uk-qm 
+vl^f {aik—\SM\bmk — q—)uk-qm 

= 6r^J^iq)fo{q) + j^tMq) - J&{q) ■ (3-48) 

Here, we define 

~ /• (fk f 1 1 A 

7 (27r)3 4^k-q U - (£^k + £^k-q) 0;+(£;k + £^k-q)y ' 

" J U-(£^k + ^k-q) ~ a; + (£^k + ^k-q)y' ' ^ ^ 
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Finally, from Eqs.(3-44) and (3-46), we obtain the following response: 

PimW = ^^imkQk ; Oab ■ (3-50) 

1 - 5^C3i7o(g) 

The pseudoscalar glueball propagator is given as P'^/e, so the pseudoscalar glueball 
mass Mq- is evaluated from the following: 

l-52c3iJo(Mo-,0) = 0. (3-51) 

Prom Eq.(3-32) and (3-51), our final task is to calculate ilo(q'), which we call the 
polarization tensor, in order to get the glueball masses. 

§4. Dependence of glueball masses on the QCD coupling constant 

The polarization tensor is rewritten in the form with the 4-momentum integra- 
tion: 

g'^no{uj,q) = no{oj,q) 

d^k 1 



(27r)3 4£;k£;k_q [U - {E^ + ^^k-q) UJ+iEk + ^k-q) 



Sk 



k"^ + ie 



The integral diverges. Thus, it is necessary to regularize the divergent integral to 
get a finite result. We here apply the dimensional regularization method. First, 
we introduce a momentum scale ^ and define the dimensionless coupling gji instead 
of g because the space-time dimension is now n. The detail calculation is given 
in Appendix A. As a result, we finally obtain the finite result for the polarization 
tensor: 

77o(,) ^ g-h,{q) = (l- ^ - -) • (4-2) 

The imaginary part of the polarization tensor leads to the decay of glueball, namely, 
it gives the decay width of glueball. However, in this treatment in this paper, only 
gluon is contained in the theory. Then, it is only possible that the glueball can decay 
to the color-octet gluon, which should be forbidden due to the color confinement. 
Thus, we here neglect the imaginary part of the polarization tensor by hand because 
the color confinement is not fully taken into account in this treatment. It is shown in 
Appendix B that the imaginary part of the polarization tensor gives a rather large 
decay width under the same parameter set. Thus, under omitting the imaginary 
part, the scalar glueball mass Mq+ and the pseudoscalar glueball mass Mq- are 
determined by the following relations given in Eqs.(3-32) and (3-51), respectively, as 

l-a±C3i7o(Mo±,0) = , (4-3) 
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1 2 3 4 5 

Fig. 1. The mass ratio between the scalar (Mq+ ) and pseudoscalar (A/q- ) gluebaU masses is shown 
as a function of the QCD running coupling qqcd ~ (?H(/i)^/(47r). The dotted lines represent 
the ratio 1 and 1.5, respectively. 



where a+ = 2 for 0^ and a_ = 1 for 0^ glueball. From Eq.(4-2) under ignoring the 
imaginary part, the glueball masses are derived as a function of the QCD coupling 
constant gn as 



Mq+ = lie ■ e 



4:77 

3a% 



A 



QCD 



exp 1 + 



lie ■ e^^R = Aqcd ■ exp 1 + 



47r2 



1 1 

3 4^ 

2 1 

3 ^ AttHo 



(4-4) 



(4-5) 



where the renormalization-group-invariant QCD scale parameter Aqqb is introduced. 
Although it may not be necessary to regard as a running coupling in the non- 
perturbative variational treatment, it is known that the running coupling can be 
derived in the lowest order approximation in this variational formalism. ^"^^ Thus, 
let the QCD coupling constant gniii) be regarded as a running coupling constant 
depending on the momentum scale fi: 



^QCD = f^e 



9R{^J) 



1 



bo In 
2 A 



^ 

QCD 



1 

8^ 



UiVc - 2iV 



/ 



.(4-6) 



Here, Nc and Nj are the number of color and flavor, respectively. We here adopt 
Nc = 2> and Nj = for the pure gauge theory, which leads to 6o = ll/Svr^. 

In Fig.l, the mass ratio is shown as a function of the QCD coupling oqcd = 
5i?,(Ai)^/(4vr). The glueball masses themselves are shown in Fig. 2 as a function of 
oqcd- Here, the QCD scale parameter is adopted as ^qcd = 0.20 GeV. For example, 
if the QCD coupling aqcD is roughly taken as oqcd = 1-6, then the glueball masses 
are obtained as 



Mo+ = 1.50 GeV 



Mn 



2.88 GeV 



(4-7) 



We cannot compare these values with experimental meson masses directly because 
the glueball must mix the scalar or pseudoscalar gg-mesons, while there exist glueball 
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3 



4 



5 



Fig. 2. The scalar glueball mass (A/q+) and the pseudoscalar glueball mass (Mq-) are shown as 
a function of the QCD running coupling Qqcd ~ 3^(^)^/(471"). The QCD scale parameter is 
taken as Aqcd = 0.20 GeV. 

candidates in particle data. So, the glueball masses obtained in our framework 
should be compared with the lattice QCD calculation such as Mq+ = 1.71 GeV and 
Mo- = 2.56 GeV.^^) Further, if we take Aqqd = 0.25 GeV, then, Mo+ = 1.53 GeV 
and Mq- = 2.58 GeV for aqcD = 2.0 are obtained. Roughly speaking, the result in 
this paper is not so bad. 



The scalar and the pseudoscalar glueball masses have been investigated in the 
framework of the time-dependent variational method and the linear response the- 
ory. We have started with the Hamiltonian of the ■sm(3) Yang-Mills gauge theory 
without fermions, namely QCD Hamiltonian without quarks. The time-dependent 
variational method within the Gaussian wavefunctional, which includes the mean 
field and quantum fluctuations around it, has been formulated in order to evaluate 
the glueball mass. The glueball mass has been calculated as the pole mass of the 
propagator of glueball. Here, the glueball propagator has been derived from the re- 
sponse with respect to the external composite field representing the glueball, which 
consists of two massless gluons. 

The gluon mass itself is zero in this method. Thus, it is shown that the finite 
glueball mass is properly generated through the interaction between massless gluons. 
Further, since the dependence of the glueball masses on the QCD coupling constant 
g reveals the form l/g"^, the results may not be arrived by the perturbation theory 
with respect to g. 

In this paper, the coupling dependence of glueball masses was given. In the 
renormalization group calculation, the glueball masses mc are expressed as 



where cq is a constant. ^'^^ Further, the similar expression was obtained in the context 



§5. Summary and concluding remarks 





18 



Y. Tsue 



of the asymptotic limit in the lattice gauge theory.^^^ In the strong coupling expan- 
sion of the lattice QCD, it was obtained that the glueball masses decrease when the 
coupling g decreases, while the change of mass ratio between O"*" state and 1"^ or 
2"'" state, instead of 0^ state, is similar to our result, namely the mass ratio increases 
when g decreases. ^^-^ Further, in the large Nc limit in the gauge/string duality, it 
seems that there is a tendency that the glueball masses decrease as g decreases. 
These results seem to be different from our result in which the glueball masses in- 
crease when g decreases. However, it may be natural that, in our result, the glueball 
masses become very large in the asymptotic region with very small g, namely in the 
quark-gluon phase, because it may be impossible that the glueballs are excited and 
produced in the deconfined phase. The investigation of the implication to the strong 
coupling expansion or gauge/string duality is an interesting future problem. 

Experimentally, it is difficult to extract the glueball masses properly because the 
glueball states mix the other (gg)-meson states with same quantum numbers. Thus, 
the glueball masses arc not fixed at present. The glueball masses obtained here 
have been compared with the results obtained by the lattice QCD simulation. The 
reasonable results are included under a certain strength of QCD coupling constant. 
However, it should be necessary to investigate other glueball states such as 2+ state. 
In addition to the glueballs with the other quantum numbers, it is interesting to 
study the excited glueball states. In order to calculate the excited glueball masses, 
the other trial states may be introduced in which {(p'\<P) = should be satisfied. 
This treatment is similar to that of the Hartree-Fock method to calculate the excited 
states in the nuclear many-body problem. Another possibility is to consider the 
three gluon states, ^^-^ where the external source term consists of three gluons. These 
investigations rest future problems. Further, in this paper, the glueball masses at 
zero temperature were considered because /„ is adopted as 1/2 in Eq.(2-17) or (2-19). 
However, if the eigenvalue of the reduced density matrix is calculated in the finite 
temperature T in which f^ = o'- (n^ + -'-/^) obtained where = l/{e'^^^/'^ — 1) 
is the bose distribution function, then, it may be possible to evaluate the glueball 
masses at finite temperature. These are future problems. 

Acknowledgements 

The author would like to express his sincere thanks to Professor K. lida. Dr. E. 
Nakano, Dr. T. Saito and Dr. K. Ishiguro whose are the members of Many-Body 
Theory Group of Kochi University. The author also would like to express his sincere 
thanks to the late Professor Dominique Vautherin for the collaboration and giving 
him the suggestion for this work developed in this paper. He is partially supported 
by the Grants-in-Aid of the Scientific Research (No. 23540311) from the Ministry of 
Education, Culture, Sports, Science and Technology in Japan. 



Scalar and Pseudoscalar Glueball Masses 



19 



Appendix A 

Evaluation of the polarization tensor nQ{uj,q) in the dimensional 
regularization scheme 



Let us show the polarization tensor in Eq.(4-1) again: 
d^k 1 



(27r)3 4E^Ek-ci 



(Ek + ^k-q) U> + {Ek + £;k_q) 



J i(2^^''^''~'' ' 



(A-l) 



Sk = 



fe2 



l€ 



Here, the 4-momentum integration is rewritten in the following form by using the 
Feynman parameter formula: 



1 



-9 



)4 P +ie(k- qY + ie 
d'^k 1 



J <27r 



)^ [{k - q^x + k'^il - x)Y 



(A-2) 



The integration of the right-hand side diverges. Thus, one needs to regularize the 
divergent integral to get the finite result. In this paper, the dimensional regulariza- 
tion method is applied and so-called modified minimal subtraction (MS) scheme is 
adopted with the consistency to the evaluation of the QCD running coupling con- 
stant Thus, the integration is regarded as the n-dimensional integration and 
is calculated as 

"1 



no{u,q) = -g 



dx 



d^k 



1 



i(2^)" [{k-qfx + k'^{l-x)Y 



g^r{2 - n/2) 
(47r)"/2 



dx- 



[q'^x{x — 1 



i2-n/2 



rl (^2)2-n/2 ^_j)2-n/2 



(A-3) 



where r{z) is the Gamma function. Here, we introduce a momentum scale ji to 
define the dimcnsionlcss coupling as follows. In n-dimension, QCD coupling g has 
a dimension. Thus, the dimensionless coupling gn should be introduced. Then, we 
can further rewrite the above result as 

,2\ /•! 



■(47r)2 



/ 



dx 



in 



x{l — x) 



(A-4) 
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where we define the dimensionless coupUng^^^ as 

e 



(A-5) 



Of course, if n = 4, then e = 0. Therefore, for infinitesimal value of e, we get 

r(e) = i-7 + 0(e) , 
e 

(AttY = i + eln47r + 0(e2) , 



(-1)-^ = 1 + ievr + 0(e2) 



2\ ^ 



= 1 - eln 



dx 



1 



- x)Y 



[ dx[l-e ln[x(l - x)]] +0{e'^) 
Jo 



= l + 2e + 0(e2) . 
Thus, the polarization tensor iTo(g) = g^no{q) can be evaluated as 



no{q) 



167r2 



- ( ^ - 7 + ln47r^ + In - ztt + 0(e) 



(A-6) 



(A-7) 



Since we apply the MS-scheme in order to get a finite value by subtracting the 
divergent term, we subtract a part proportional to the following set: 



- = 7 + In 47r . 

e e 

Finally, as a result, we get the polarization tensor as 

.2 ^2 



Appendix B 

Decay width 



(A-9) 



In this Appendix, the glucball mass is reconsidered by taking into account the 
imaginary part of the polarization tensor which may lead to the decay width of 
the glueball. In general, when a mass function U{p) has an imaginary part, the 
propagator has a form 

-1 



Sij cx 



P 



p2 - M2 - 



+ iMr 



(B-l) 



where the mass function E is divided into a real and an imaginary part as U = 
M - ir/2. Thus, the relations 



Re ^"^(0;,^ = 0) = , 
Im 5-\a;,p = 0) = MF 



(B-2) 
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1 2 3 4 5 " 

Fig. 3. The scalar glueball mass (Mq+) and the pseudoscalar glueball mass (Mq-) are shown as a 
function of the QCD running coupling oqcd ~ <?fl(M)^/(47r), in which the imaginary part of the 
response function is taken into account. The QCD scale parameter is taken as Aqcd ~ 0.20 
GeV. 



give the mass M and F that may be regarded as a decay width. 

In the treatment developed in this paper, the above relations are rewritten for 
the scalar 0+ and pseudoscalar 0^ glueballs as 



1 — a^csRe 77o(wo±, 0) = , 

Mo±ro± 



^0± 



M, 



o± 



r2 



-a±C3lm no{ujQ±,0) 



0± 



with a+ = 2 for 0+ and a_ = 1 for glueball. Here, 

^2 



167r2 



ITT 



From Eq.(B-2), M and F are obtained as 



(B-3) 



(B-4) 



fie ■ e^^i^'fl 



^QCD • exp (l + ^ (glj + ^ 



\^ Wit J 
aqCD (j 



3a±«QCD 



-j^ _|_ ^ 3a±aQCD 



8Mr 



0± 



3a±aQCD 



1^ _^ ^3a±aQCDy 



(B-5) 



(B-6) 
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In Fig.3, the glueball masses are shown as a function of oqcd = Qr/'^t^ with vIqcd = 
0.2 GeV. If the QCD coupHng constant oqcd is roughly taken as agcD = 2.0, then 
the glueball masses are obtained as 

Mo+ = 1.76 GeV , Mq- = 2.44 GeV . (B-7) 

Under the above parameter set, from the imaginary part of the polarization tensor, 
the decay width may be evaluated, which results r'o+ = 2.54 GeV and Fq- = 2.61 
GeV. These values are rather large, while a large decay width is reported by using 
a chiral quark model. As is mentioned in §4, only gluons are contained in this 
framework. Then, it is only possible that the glueball decays to color-octet gluons, 
which should be forbidden due to the color confinement. However, in our treatment, 
the color confinement is not considered explicitly, so the glueballs easily decay to 
gluons. Thus, the rather large decay width, which means a decay from the color- 
singlet glueball to the color-octet gluons, may be obtained unavoidably. Thus, the 
investigation of the decay of glueball is still remained as an open question. 

Appendix C 

Gluon mass 

In this appendix, we show that the gluon mass itself is zero in the framework 
developed in this paper under the Gaussian approximation used there. First, the 
canonical equations of motion for yl"(a;,t) and E'^{x,t) in Eq.(2-9a) with (2-7) are 
given as 

= -e,jkdM{x) + ge,jkr'^B]{x)Al{x)-\ J d3^^^GFj'(y,y), (C-1) 
where we define 

geiimemjkr" {-Sdbdl - gf^'^Aliy, t)) 6\x - y) 

gejimemikf'"' {-Sdadl - gf^'Aliy, t)) S^{x - y) . (C-2) 

In order to get the gluon propagator, we have to introduce the external term i^ext 
in the Hamiltonian as 

H^^t = J d^xJ^{x,t)A1{x) = J^e-''^^ j (fxe^'i-'^A1{x) . (C-3) 

Thus, the gluon propagator 5?^ can be derived by 6{A^)/6Jj following the general 
discussion. 



SAi{x,t) 
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Let us start with solutions A'^(x,t) = El{x,t) = under the Hamiltonian Hq. 
With the external term, the solutions should be shifted. Here, we denote them as 



SA^ix) 



(C-4) 



From the equation of motion in (C-1), we can get the equation of motion for dA'^(x, t) 
witli a Unear approximation for dA^ {x,t) under small source current J": 

{{d^ - V^)Sij + didj) SA';{x,t) - ^g\daSjk - 26ikdji + 5ijdM)Gij{x)SAl{x,t) 



(C-5) 



where we introduced a new notation Gij{x) through Gfj{x,x) = SabGij{x). In the 
above equation of motion, the second term, Gij{x), is diagrammatically represented 
by so-called tadpole diagram. It is well known that there is no tadpole contribution 
in pure Yang-Mills gauge theory in the dimensional regularization scheme, namely. 



Gij{x) oc 



f d^k 


1 


J (27r)3 2|fc| 


/• d''k 




J «(27r)" 




1 r 




(47r)t 1- 


n 
2 


1 




(47r)2 


[a 



d'^k 



-1 



i(27r)4 kl - |fe|2 + ie 



m2-)-0 

7 1 — In -|- 0(e) 







(C-6) 



Thus, from the equation of motion in Eq.(C-5), the following equation in the mo- 
mentum space is obtained: 



iJ'Sij -\q\ 



l«|2 



(C-7) 



Finally, the gluon mass is given by the pole of the gluon propagator Sfj{oj,q) = 
SA'^/Jj, namely, gluon mass is exactly zero in this framework. 
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